Ferromagnetism in the Two-Dimensional Periodic Anderson Model 
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Using the constrained-path Monte Carlo method, we studied the magnetic properties of the two- 
dimensional periodic Anderson model for electron fillings between 1/4 and 1/2. We also derived 
two effective low energy theories to assist in interpreting the numerical results. For 1/4 filling we 
found that the system can be a Mott or a charge transfer insulator, depending on the relative values 
of the Coulomb interaction and the charge transfer gap between the two non-interacting bands. 
The insulator may be a paramagnet or antiferromagnet. We concentrated on the effect of electron 
doping on these insulating phases. Upon doping we obtained a partially saturated ferromagnetic 
phase for low concentrations of conduction electrons. If the system were a charge transfer insulator, 
we would find that the ferromagnetism is induced by the well-known RKKY interaction. However, 
we found a novel correlated hopping mechanism inducing the ferromagnetism in the region where 
the non-doped system is a Mott insulator. Our regions of ferromagnetism spanned a much smaller 
doping range than suggested by recent slave boson and dynamical mean field theory calculations, 
but they were consistent with that obtained by density matrix renormalization group calculations 
of the one-dimensional periodic Anderson model. 
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I. INTRODUCTION 

Identifying the origin of itinerant ferromagnetism in 
metals and specifying simple models exhibiting it are two 
of the most intriguing and long-standing problems in con- 
densed matter physics. Here we report the results of low- 
energy perturbation theory calculations and supporting 
zero temperature quantum Monte Carlo (QMC) simula- 
tions that suggest the existence and mechanisms for fer- 
romagnetic (FM) ground states in the two-dimensional 
periodic Anderson model. 

/^From a historical point of view, the one-band, nearest- 
neighbor hopping Hubbard model was one of the first 
models proposed to describe itinerant ferromagnetism; 
however, the ferromagnetic phase has never been found 
at physical parameter values. The numerical calcula- 
tions, for example, have narrowed the extent of this phase 
to a small region around the Nagaoka point pi, that 
is, the strong coupling limit for one hole doped away 
from half-filling. Paradoxically, in two-dimensions (2D), 
this model exhibits anti-ferromagnetism at half-filling 
and anti-ferromagnetic correlations around half-filling at 
weak and intermediate couplings. 

Recently, Guerrero and Noack B listed several possi- 
ble extensions of the Hubbard model that should enhance 
ferromagnetism: (i) the addition of frustrating hopping 
terms [&-p|, (ii) the inclusion of more than one orbital per 
unit cell, and (iii) the addition of more general nearest- 
neighbor interactions [p|q|. In fact, a number of frus- 
trated models with more general interactions, such as 
the t — t' Hubbard model, and multiband models, such 



as the periodic Anderson model (PAM), have ferromag- 
netic ground states pBM]. In this paper, we focus on the 
properties of the two-dimensional PAM. 

The PAM is often used to describe the essential 
physics of many transition metals and rare-earth and ac- 
tinide metallic compounds including the so-called heavy- 
fermion systems g . The model includes a band of "light" 
uncorrelated electrons hybridized with a band of heavy 
strongly correlated electrons. Despite intense efforts to 
determine its phases, only a few controlled analytical ap- 
proximations and numerical calculations exist for D > 2. 
Instead, previous work often studied the single-impurity 
Anderson model and focused on the competition between 
Kondo screening and the direct RKKY coupling between 
the localized spins |^,0- This competition is present 
when the number of conduction electrons is at least sim- 
ilar to the number of singly-occupied low-lying /-levels. 
In D > 2, much less attention has been given to the re- 
gion of the phase diagram where the density of electrons 
in the conduction band is small; however, it is for this 
case that ferromagnetism has been established in one di- 
mensional systems |7|JllHl5| . 

Using the density matrix renormalization (DMRG) 
method in the one-dimension, Noack and Guerrero [fj, 
for example, found partially and completely saturated 
ferromagnetism in the PAM. They considered a param- 
eter regime where the energy of orbital state £f and the 
strength of the Coulomb repulsion U were adjusted so 
each orbital had just one electron. The position of the 
orbital energy was below the lower band of the non- 
interacting problem, and one electron per orbital cor- 



responds to a 1/4-filled non-interacting problem. For a 
sufRciently large value of U, the model exhibited a ferro- 
magnetic ground state. Beyond an interaction-dependent 
value of the doping and a doping-dependent value of [/, 
this state disappeared. The ferromagnetic phase was a 
peninsula in a phase diagram that was otherwise a sea of 
paramagnetism except at 1/4 and 1/2 filling where the 
ground state of the PAM was antiferromagnetic. 

Ferromagnetism seems to be readily found by mean- 
filed approximations in any dimensions |pq-E2|. Using a 
slave-boson mean-field theory (SBMFT) for the symmet- 
ric PAM, Moller and Wolfe |ll6| found results similar to 
those of Noack and Guerrero. At 1/2 filling, they found 
a paramagnetic (PM) or antiferromagnetic (AF) phase 
depending on the value of the Coulomb repulsion U. By 
lowering the density of electrons from 1/2 filling, they 
also found a smooth crossover from AF to FM order via 
a spiral phase. Just before 1/4 filling, they got a first- 
order transition from FM to AF order. More recently, the 
SBMFT calcualtions of Doradzihski and Spalek [0jl8| 
found wide regions of ferromagnetism in the intermedi- 
ate valence regime that surprisingly extended well below 
1/4 fiUing. 

In the low temperature dynamical mean-field theory 
(DMFT) calculations, Tahvildar-Zadeh et al. ^ also 
found a region of ferromagnetism and studied its temper- 
ature dependence. At very low temperatures, their ferro- 
magnetic region extended over a wide range of electron 
density and in many cases embraced the electron density 
of 3/8. At 3/8 filling they proposed a specific Kondo- 
induccd mechanism for ferromagnetism that has the con- 
duction electrons in a spin-polarized charge density-wave 
anti-aligned with the ferromagnetically aligned local mo- 
ments on the valence orbitals. More recently Meyer and 
Nohing l^-ll appended perturbation theory to DMFT 
and also predicted ferromagnetism over a broad ramge of 
electron filling extending below 1/4 filling. 

Our previous [g3[ and new QMC results qualitatively 
agree with the DMRG work; however, the phases we 
find quantitatively and qualitatively disagree with those 
derived from the mean-field approximations. Quanti- 
tatively, we find ferromagnetism in a narrower doping 
range than the one predicted by the DMFT and SBMFT 
calculations. For densities between 3/8 and 1/2, QMC 
predicts a paramagnetic region, whereas mean-field the- 
ory predicts a ferromagnetic states in part of that re- 
gion. In fact, at a filling of 3/8 where these calcula- 
tions predict ferromagnetism, we find a novel ground 
state of an entirely different symmetry. Instead of fer- 
romagnetism, QMC finds a resonating spin density-wave 
(RSDW) state; that is, the ground state was a linear 
combination of two degenerate spin-density waves char- 
acterized by the (tt, 0) and (0, n) wave vectors. 

We remark that the quantitative differences between 
the DMRG and QMC calculations and the DMFT and 
SBMFT calculations most likely result from the expected 



breakdown of mean-field theory in one and two dimen- 
sions. Probably the RSDW state was not found because 
it was not sought. On the other hand, trying to under- 
stand the mechanism for ferromagnetism is more funda- 
mental. It points to the long-standing difficulty of build- 
ing an understanding of the PAM upon the better un- 
derstood single impurity Anderson model (SIAM) or the 
analogous problem of building an understanding of the 
Kondo lattice model upon the better understood single 
Kondo impurity problem. Nozieres addressed this later 
problem and proposed a reconciliation in terms of what 
he calls "the exhaustion picture" |2^,|2^ . 

We note that the electron densities near 1/4 filling 
place the work of Noack and Guerrero (t) in the exhaus- 
tion regime. When the Coulomb repulsion U associated 
with the double occupancy of an orbital is large (strong 
coupling), it is often argued that around 1/4 filling the 
PAM behaves like a 1/2-filled one-band Hubbard model. 
In two-dimensions this Hubbard model has an antiferro- 
magnetic ground-state generated by a superexchange in- 
teraction Jsx ^ V^t'^/U^ (for the symmetric case) where 
V is the strength of the hybridization in the PAM. It 
has also been established that the Kondo lattice model 
and the Hubbard model for strong coupling are isomor- 
phic with particles in one becoming holes in the other 
pTI . The relevance of the Hubbard model gains addi- 
tional support from Nozieres's long-standing argument 
that the dynamics of the screening clouds are described 
by an effective Hubbard model. 

Because of these mappings, it is seems consistent 
to suggest the strong-coupling physics found by Noack 
and Guerrero would map onto the domain of a one- 
dimensional 1/2-filled Hubbard model. The venera- 
ble theorem of Lieb and Mattis ||2^, however, excludes 
the possibility of the Hubbard model (with nearest 
neighbor hopping) from showing ferromagnetism in one- 
dimension. Their proof relied on the obvious ability to 
order the electrons along the chain. The two-band nature 
of the PAM, however, prevents this ordering by allowing 
processes not possible in the Hubbard model. A similar 
situation would occur for the two-legged Hubbard model, 
if it in fact shows ferromagnetism. We will argue that in 
the two-dimensional PAM these same processes are re- 
sponsible for the ferromagnetism. They are RKKY and 
Nagaoka-like and are excluded in Nozieres's [^,g5| pic- 
ture and Tahvildar-Zadeh et a/.'s [ |l9[ interpretation of 
it. It is important to remark that ferromagnetism in the 
PAM is obtained for a large region of parameters which 
include realistic values. 

We will base our arguments on the predictions of effec- 
tive Hamiltonians generated from the PAM by perturba- 
tion theory and the results of our QMC simulations. For 
the PAM parameters studied, these effective Hamiltoni- 
ans suggest a paramagnetic or anti-ferromagnetic state 
at 1/4 filling, a ferromagnetic region between 1/4 and 
3/8 filling, a RSDW at 3/8 filling, a paramagnetic region 



between 3/8 and 1/2 filling, and an antiferromagnetic 
state at 1/2 filling. We see all these features in the QMC 
simulations. 

In the Section II we will define the PAM and sketch our 
derivation of the effective Hamiltonians. In Section III we 
will summarize our numerical method, noting finite size 
limitations. Our results will be presented in Section IV. 
In the Section V, the Concluding Section, we will give 
a detailed contrast between our picture and select other 
works. 
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The PAM is described by the Hamiltonian 
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where dj^ and //^ create an electron with spin a in d and 
/ orbitals at site i in a square lattice, and nj^ — fl^fi^ 
is the number operator for the /-electrons of spin a at 
site i. Elsewhere we will use a similar notation to de- 
note quantities like nf^ = d\^d^^, the number operator of 
(i-electrons. The lattice has N sites and the hopping am- 
plitude t between d-orbitals is only to nearest-neighbor 
(n. n.) sites. The hopping amplitude V hybridizes differ- 
ent orbitals on the same site. We used periodic boundary 
conditions. 

^From (nl) we define Hq, the resulting Hamiltonian 
when [/ = 0. Hq has two dispersive bands 
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separated by a gap 
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FIG. 1. Band structure of the non-interacting (f/ = 0) 
two-dimensional periodic Anderson model. In units of f, 
V = -0.5. In (a) e/ = 0.2 and (b) e/ = 0.5 Ilustrated is 
the flatnes of the lower band for case (b). 

fS) ^^ ^^^° note that the widths of the upper and lower 
band are 



For a square lattice, the energy Ck ~ — 2t(cos k^+cos ky). 
This band structure for Hq is illustrated in Fig. n^. We 
note if £/ becomes very negative (doping way below the 
bottom of the lower band), A approaches |e/|. 
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As 6/ becomes very negative, W approaches zero and 



W~^ approaches W — 8t, the band width when U = V = 
0. 

The operators which create quasi-particles in the lower 
and upper bands are of the form 



^ka 



with 



Uk 



Vk 



^k/L + «k4^ 

/^L = -%/L + "k4<T 



(5) 
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The symmetric PAM, which has the electron filling 
p = 1/2 and U = — e/2, has particle-hole symmetry. 
This symmetry is sufficient to prevent the fermion sign 
problem that plagues QMC simulations. Such simula- 
tions, performed by Vekic et al. [g8[, suggest the exis- 
tence of a charge and spin gap for small values of U with 
the spin gap disappearing when U is increased to some 
Uc ~ 2. Above Uc the system exhibits long-range anti- 
ferromagnetic order. 

In a previous work we presented QMC results for the 
asymmetric model. QMC simulations for the asymmetric 
model experience a sign problem which is the reason why 
we used the constrained-path Monte Carlo method [ p9[ . 
For fixed values of e/ , we varied U and hole-doped away 
from half-filling down to a filling of 3/8. For a large 
enough value of U we also found anti-ferromagnetism 
at 1/2 filling. This state was rapidly suppressed upon 
doping. At 3/8 filling we saw a sharp peak appear in 
the spin-spin correlation at the wavenumbers k = (0, tt) 
and (0, tt). We interpreted this peak as consequence of 
a state resonating between two degenerate spin-density 
waves characterized by reciprocal wavevectors (0,7r) and 
(0,7r). 

In the present work we explore the doping range from 
1/4 to 3/8 filling, arguing for a region of ferromagnetism. 
Part of our arguments will be based on the properties 
of effective Hamiltonians for two different regions of pa- 
rameters used in the simulations. These Hamiltonians 
will be derived in Section IV. We now summarize the the 
constrained-path Monte Carlo method. 



III. NUMERICAL METHOD 

The constrained path Monte Carlo (CPMC) method 
is extensively described and benchmarked elsewhere [p9| . 
Here we only discuss its basic strategy and approxima- 
tion. In the CPMC method, the ground-state wave func- 
tion lipo) is projected from a known initial wave function 
IV't) by a branching random walk in an over-complete 
space of Slater determinants |0). In such a space, we can 



write IV'o) = J2(t> x{^)\4')- The random walk produces an 
ensemble of {(p), called random walkers, which represent 
\ipo) in the sense that their distribution is a Monte Carlo 
sampling of x((/>), that is, a sampling of the ground-state 
wave function. 

More specifically, starting with some trial state \iPt), 
we project out the ground state by iterating 



m = e 



_ -At{H-Et] 



W 



(7) 



where Et is some guess of the ground-state energy. Pur- 
posely At is a small parameter so ioi H = T + V we can 
write 



e-^^" « e 



-AtT/2 -AtV -AtT/2 



(8) 



where T and V are the kinetic and potential energies. 

For the study at hand, the initial state \iPt) is the di- 
rect product of two spin Slater determinants, i.e.. 



i^T)=ni<^? 



(9) 



Because the kinetic energy is a quadratic form in the cre- 
ation and destruction operators for each spin, the action 
of its exponential on the trial state is simply to transform 
one direct product of Slater determinants into another. 
While the potential energy is not a quadratic form in 
the creation and destruction operators, its exponential 
is replaced by sum of exponentials of such forms via the 
discrete Hubbard-Stratonovich transformation 
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provided U > and coshArJ = e" 
we re-express the iteration step as 



-AtU/2 



Accordingly 
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where x = {xi,X2, ■ ■ ■ ,Xn) is the set of Hubbard- 
Stratonovich fields (one for each lattice site), N is the 
number of lattice sites, P{x) = (|)^ is the probability 
distribution for these fields, and B„{x) is an operator 
function of these fields formed from the product of the 
exponentials of the kinetic and potential energies. 

The Monte Carlo method is used to perform the multi- 
dimensional integration over the Hubbard-Stratonovich 
fields. It does so by generating a set of random walkers 
initialized by replicating j-^x) many times. Each walker 
is then propagated independently by sampling a x from 
P{x) and propagating it with B{x). After the propaga- 
tion has "equilibrated," the sum over the walkers pro- 
vides an estimate of the ground-state wave function jV'o)- 



In practice we performed an importance-sampled ran- 
dom walk by using the transformed iterative equation 
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dxP{x)Y[B,{x)\Y[> 



In this equation 



P{x) = ZP{x) 



Y{ArT\Ba{x)\<l>.) 

n.(0Ti'/'-> 



(11) 



(12) 



Thus importance sampling changes the probability dis- 
tribution of the Hubbard-Stratonovich fields, biasing it 
towards the generation of states with large overlap with 
the initial state. The factor Z is the normalization con- 
stant for the new distribution. It is associated with the 
weight assigned to each walker and, the weight is used in 
a branching process to control the variance of the results. 
We will not discuss this process here. 

We used two different estimators for the expectation 
values of some observable O. One is the mixed estimator 



{O). 



(V'tIV'o) 



and the other is the back-propagated estimator 
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where \tpo) is the QMC estimate of the ground state and 
i is typically in the range of 20 to 40. For observables 
that commute with the Hamiltonian, the mixed estima- 
tor is a very accurate one and converges to the exact 
answer as |V'o) converges to exact ground state. For ob- 
servables that do not commute with the Hamiltonian, 
like correlation functions, the back-propagated estima- 
tor has been found to give very accurate estimates of 
ground-state properties. Significant differences between 
the predictions of these two estimators often exist. 

To completely specify the ground-state wave function 
for a system of interacting electrons, only determinants 
satisfying {iJo\4'a) > are needed because It/iq) resides in 
either of two degenerate halves of the Slater determinant 
space, separated by a nodal surface N that is defined by 
{^o\(j)a) = 0. The degeneracy is a consequence of both 
jf/'o) and — l^o) satisfying Schrodinger's equation. The 
sign problem occurs because walkers can cross N as their 
orbitals evolve continuously in the random walk. Asymp- 
totically they populate the two halves equally, leading to 
an ensemble that has zero overlap with IV'o)- If N were 
known, we would simply constrain the random walk to 
one half of the space and obtain an exact solution of 
Schrodinger's equation. In the constrained-path QMC 
method, without a priori knowledge of N, we use a trial 
wave function ItJjt) and require (t/itI^o-) > 0. This is 
what is called the constrained-path approximation. 



The quality of the calculation clearly depends on the 
trial wave function IV't)- Since the constraint only in- 
volves the overall sign of its overlap with any determinant 
{(p), it seems reasonable to expect the results to show 
some insensitivity to \tpT)- Through extensive bench- 
marking on the Hubbard model, it has been found that 
simple choices of this function can give very good results 

Besides as a starting point and as a condition con- 
straining a random walker, we also use |V't) as an impor- 
tance function. To reduce variance, we use ('0t|'/'o-) to 
bias the random walk into those parts of Slater determi- 
nant space that have a large overlap with the trial state. 
For all three uses of \iPt), it clearly is advantageous to 
have IV't) approximate ji/io) as closely as possible. Only 
in the constraining of the path does JV't) t^ IV'o) generate 
an approximation. 

We constructed |V't) ~ Yia I'^t) from the eigenstates of 
the non-interacting problem. Because the z-component 
of the total spin angular momentum is a good quan- 
tum number, we could choose unequal numbers of up 
and down electrons to produce trial states and hence 
ground states with Sz = ^{N^ ~ ^i)- Whenever pos- 
sible, we would simulate closed shells of up and down 
electrons, as such cases usually provided energy estimates 
with the least statistical error, but because we wanted to 
study the ground state energy as a function of 5"^, we 
frequently had to settle for just the up or down shell be- 
ing closed. In some cases, the desired value of Sz could 
not be generated from either shell being closed. Also 
we would select the non-interacting states so \iPt) would 
be translationally invariant, even if these states used did 
not all come from the Fermi sea. The use of unrestricted 
Hartree-Fock eigenstates to generate 10^ instead of the 
non-interacting eigenstates generally produced no signif- 
icant improvement in the results. 



IV. RESULTS 



Effective Hamiltonian: Wannier Orbital 
Approach 



Our first effective Hamiltonian explicitly targets cases 
where the lower band of the non-interacting model is very 
fiat. Such cases exist for — e/ > W/2 > V > 0. In this 
regime of parameters and around 1/4 filling, single elec- 
tron occupancy of the f-states can occur because of the 
depth of the orbital state as opposed to the double occu- 
pancy penalty of the Coulomb repulsion. We will begin 
by building Wannier orbital operators [30[ for each band 
from the quasi-particle operators defined in (H) 



where Ry = r^ — Tj . 
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Rewriting Hq in the Wannicr basis, we find that 



(17) 



^,j,c 



with 



{l.n) I 

+ y'^{b%iaji +bjiaii) 
I 

- V^^(6i;aj7 +bjiaii) 



(18) 



i/o is simply the sum of two hopping terms correspond- 
ing to the lower and the upper bands. By construction, 
no hybridization exists between these two bands. The 
cost for this simplification is non-zero hoppings t", and 
Tj'- between any pair of Wannier orbitals i and j in the 
same band. 

Next we rewrite the interaction term 



Hr, 



in the Wannier basis 



^-E 






(19) 



J,o- 



Hij = U ^ (fly all + fey/3||)(a,,j-a,,| + b^.^fi^,^) 

X {aijaj^ + bi^pl^){ai,jai,^ + bi,j(ii,{) (20) 

This expression appears more complex than the one in 
the original basis; however, from it we can more con- 
veniently derive a low energy effective Hamiltonian for 
electron fillings less than 1/2 filling. 

To do this we first require that A > C/, i.e. the sys- 
tem is a Mott insulator for 1/4 filling, so we can initially 
consider a H\ ^^^^^ that does not have processes involving 
the upper band 

(21) 



with gaiii' = J2j ^ij^i'j'^ij^i'j- To identify more easily 
the physically different contributions, we rewrite the U 
term to produce 
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i,i' ,1,1' 
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Jti=2UY,a. 
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a^jUjiajv 



(23) 



where ^ means that there are no repeated indices. We 
see that HI ^j^^^ is an extended Hubbard model with long 
range hoppings, a ferromagnetic exchange interaction, 
correlated hoppings, and a term which destroys a spin- 
anti-aligned pair of electrons in sites i' and I' and creates 
an anti-aligned pair at i and I. 

Again for a large range of parameters, the lower band 
of the PAM is quite flat. If we regard |e/| as very large, 
we can Taylor series expand Eq. o, substitute the result 
into Eq. Hq, and obtain 



-tV^I 



k/l' 



for i = j 

for i and j n. n.) 



-y/|e/|, fori = j 

— ty/|e/| , for i and j n. n. 



(24) 



Matrix elements for i and j beyond nearest neighbors 
(n. n.) are smaller by higher powers of y/|e/|. Thus 
the Wannier operator a- is predominately / as the am- 
plitudes aij and bij strongly decrease with the distance 
between i and j. 

With these results we see that J^^ and uouii are 
poportional to t^V^ /\ef\^, while gumi is proportional to 
t^V^ /\€f\^ . This means that we can neglect the last term 



i,i' ,1,1' 



of H 
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In this one band Hamiltonian the Coulomb repulsion 
is no longer just on-site. Its spatial extension depends on 
the spatial extension of the Wannier orbitals. Because of 
original form of this interaction, it still only affects sin- 
glet states, but these states can now be non-local. We 
note that even if U /t ^ 1 this one-band model is in the 
strong-couphng regime: tj ^ U and r^" « ^(^^/k/P) 
if i and j are nearest neighbors so U /t^, ^1. In the 
new basis the narrow band appears narrower, the heavy 
fermions appear heavier, but the interaction experiences 
little renormalization. 

The non-locality is the origin of the ferromagnetic 
Heisenberg term. In this term, provided U is not ar- 
bitrarily large, the ferromagnetic exchange interaction 
Jfi ~ Ut^V'^ /\eff' , however, is smaller than the antifer- 
romagnetic super-exchange interaction which is of order 
(Tj")^/C/ ^ t^V'^ /U\ef\'^ . we remark is that the magni- 
tude of the AF interaction is very small. For this reason, 
the lowest order terms involving the upper band are cru- 
cial to determine the magnetic phase of the system doped 
above 1/4 filling. 

These lowest order upper band processes come from 
terms in the Eq. Effl with one (3 operator and can be writ- 
ten as 



<L.-c^Ev 



(«L/3,* + /3la.ff) (26) 



with 
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(27) 



Here the terms with I ^ V can be neglected when the 
lower band is flat so we can rewrite Eq. E^ as 



hja <^ |--T(*»,i + T—M-]\.l) 



|e/l 



|e/l 



(28) 



Thus the lowest order inter-band process are correlated 
hoppings between ai and [3j orbitals and are proportional 
to the spin polarization (opposite to the spin of the elec- 
tron which hops) surrounding the j site: the hoping oc- 
curs out of regions of ferromagnetically aligned electron 
spins. 

It is interesting to note that the origin of this ferromag- 
netic alignment is not an exchange mechanism but pro- 
cesses involving charge-transfer. To see this more clearly 
we show in the Appendix that by using a Schrieffer- Wolff 



transformation we can reduce the multi-band term to an 
effective one band term 
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with Sk ~ cos fcj. + cos ky . The hopping between two a 
orbitals is only possible if both sites are occupied with 
electrons having the same polarization. Then it is clear 
that H\^ hand induces a ferromagnetic interaction between 
the localized /-state which comes from the itinerancy of 
the carriers trough the upper band. In addition, to max- 
imize the energy gain, the added carriers must occupy 
k ^ states. In this way, the magnitude of the hopping 
tij(k) is maximized at the same time the denominator 
_E+ (k) — tf — Uis minimized (tends to A — [/) . The elec- 
trons can be added to k ^ states only if the background 
is ferromagnetic and, of course, opposite to the spin of 
the added electron. These charge-transfer processes en- 
hancing ferromagnetism involve the states in the lower 
part of the upper band. 

This inter-band process is illustrated in Fig. ya. The 
representation emphasizes the collective nature of the 
mechanism. The virtually hopping electron has reached 
a k = band state. It is anti-aligned with the /-orbitals 
which are ferromagnetically aligned among themselves. 
The moment of this band state partially compensates the 
fully saturated ferromagnetic alignment of these orbitals. 
This compensation picture differs from "the exhaustion 
picture" of Nozieres 24,23] evoked by Tahvildar-Zadeh et 
al. |M The compensation is on a collective state to col- 
lective state basis and not the collective state to single 
moment basis argued by Nozieres. This difference high- 
lights the difffculty building the physics of the periodic 
Anderson model from the physics of its single impurity 
version. 

The process in Fig. ||b contrasts that in Fig. ga. This 
process in Fig. gb compensates on a site-to-site basis 
and is the one present for a paramagnetic or an anti- 
ferromagnetic alignment of the /-orbitals. The energy 
cost for virtual hopping is higher in than that in Fig. Ha. 
This leads to larger energy denominators in perturba- 
tion theory and in turn leads to a smaller lowering of the 
ground state energy. 
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b) Non-magnetic Background 

FIG. 2. Conduction electron compensation of /-orbital mo- 
ments, (a) The mechanism for the effective one-band Hamil- 
tonian. An electron at the bottom of the conduction band 
is partially compensating the coUecive ferromagntic state. 
This process more effectively lowers conduction electron ki- 
netic energy than mechanism (b). In the latter the conpen- 
sation is a one-to-one on-site process and is present for an 
anti-ferromagneitc or paramagnetic aligment of the /-orbitals. 
A paramagnetic alignment is shown 



B. Effective Hamiltonian: Canonical Transformation 

Here we will present an effective low energy model 
valid in a different region of parameter space: U 3> 
14^/2, —£f ^ V > 0. In this regime of parameters and 
around 1/4 filling, single electron occupancy of the /- 
states occurs mainly because of the double occupancy 
penalty of the Coulomb repulsion. As we will see, the in- 
teraction between the moments in these /-states is dom- 
inated by the RKKY interaction. 

To derive this effective Hamiltonian, we will make a 
fourth order Schrieffer- Wolff transformation [Q as in 
Ref. M: 



H^e^'e^'He' 



W1-W3 



(31) 



where the transformation operators Wi and W3 are of 
order V and V^ . With this transformation we get a new 
Hamiltonian H without terms of order V and V^. By 
means of another canonical transformation we eliminate 
the term of order V^ . In this way we get the low energy 
effective Hamiltonian Hspm correct through V*. The de- 
tails of the derivation are given in the Appendix. 



The final expression for the effective Hamiltonian is a 
Heisenberg Hamiltonian 
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The RKKY contribution J^- ' is given by 
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)o is defined in the Appendix. 

(1) 



Jj is associated with virtual processes where one /- 
electron at sites i and j go through the conduction band, 
doubly occupy the d-sites j and i, and one of these two 
electrons comes back to an empty /-orbital. The expres- 
sion for J;, is 
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The other contributions to Jij are given by 
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where m-kk', f^kk', 2:kk' and 2/kk' are defined in the Ap- 
pendix. These four contributions to Jy are about one 
order of magnitude smaller than J^- . 

We expressed the magnitude of the effective spin-spin 



interactions Jij ( |33| ) in terms of k-space states because 
in the absence of hybridization they are the eigenstates 
of the problem. From this point of view, the ferromag- 
netism for low density of carriers comes from the RKKY 
term J^ and the small volume enclosed by the Fermi 
surface: When the Fermi volume is small, the transferred 
k-vector between different points on the surface is small, 
and the relevant phase factors of Eq. p3 are positive for 
small distances between sites i and j. This positivity 
gives an effective ferromagnetic interaction between near 
neighbors, which is proportional to the density of carriers 



and competes with the usual anti-ferromagnetic super- 
exchange interaction included in the other components of 
Jij . Because the ferromagnetic interaction increases with 
the concentration of carriers, we expect the appearance of 
a ferromagnetic phase above some critical concentration 
Pc- Due to the non-interacting case {U = 0) being para- 
magnetic, we also expect a ferromagnetic phase above 
some critical value Uc- To derive a simple expression 
for the phase boundary Uc = U{pc), we have to com- 
pare the effective ferromagnetic interaction J^ ■ with the 
the antiferromagnetic super-exchange interaction, which 
is proportional to V'^t'^ /U^ for the symmetric PAM. From 
Eq. ^ it is clear that, for the dilute case, j|- is propor- 
tional to the concentration of carriers pc and V'^/tU'^. 
Consequently, the phase boundary of the ferromagnetic 
region is given by 



PcV^/tu^ ~ vHyu'' 



which implies that 



Pc 



tyu! 



(40) 



(41) 



This expression qualitatively agrees with the numerical 
results of Ref. @. 

This ferromagnetic phase must disappear above some 
upper critical concentration of carriers where the vol- 
ume enclosed by the Fermi surface is no longer small. 
For this concentration we do not expect an ordered 
phase unless nesting is present. In presence of nesting, 
_B~(k) = E~{]s.+ Q) = i?Formi, the Fermi surface may 
be unstable towards the development of a spin density 
wave. This instability is manifested in Eq. |34| from the 
divergence of the Q Fourier component of J^ ■ . For a 
commensurate state, Q is a high symmetry point of the 
Brillouin zone and equals one-half of a reciprocal lattice 
vector. For a square lattice two such points exist: (0,7r) 
for p = 1/4 and (tt, tt) for p = 1/2. In the absence of nest- 
ing no well defined Q vector connects different points of 
the Fermi surface, and preferred ordering for the spins is 
absent. For this reason, if we were to increase p over the 
upper ferromagnetic phase boundary, we would expect a 
paramagnetic phase for any concentration different from 
1/4 and 1/2. This picture is similar to the one already 
obtained in ID [Q. The only difference between ID and 
2D is in ID nesting exists at the Fermi level for any con- 
centration of carriers. This nesting explains why spiral 
critical order at 2\ip is obtained for any value of p in ID 



and 12 X 12 one-band Hubbard model. As mentioned 
previously, we used ground states derived from the non- 
interacting problem as \iJjt)- When we simulated the 
PAM with £f = —5, we used the non-interacting states 
for 6/ = — 2 because this choice consistently produced a 
lower estimate of the ground state energy with a smaller 
statistical error then we would obtain if we had used the 
non-interacting states for e = —5. 

Most of our simulations were performed for electron 
fillings of 1/4 through 3/8. In our previous work P3| , 
where we studied the region from 3/8 to 1/2 and found 
the RSDW ground state at 3/8 filling for a 6 x 6 sys- 
tem. At 1/2 filling, we found the expected AF ground 
for [/ > 2 in agreement with other QMC simulations. In 
between, we found a paramagnetic (PM) ground state. 
For these previous studies, we had e/ = — 2. As part 
of the present study, we repeated some simulations over 
this 3/8 to 1/2 range with e/ = —5 but did not find any 
indications of a FM state. We note again that the predic- 
tions of DMFT and SBMFT are inconsistent with these 
PM and RSDW states. The effective Hamiltonians pre- 
sented in this paper are expected to be inappropriate for 
this filling range. We will not discuss this range further. 

In the present work, we also did a series of rough sim- 
ulations at fillings less than 1/4 for various values of U 
and £/. We only found PM ground states. DMFT and 
SBMFT find FM ground states for identical parameter 
choices. 
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FIG. 3. Energy difference between the partially polarized 
and singlet ground state energies as a function of electron fill- 
ing. The lattice has 6x6 unit cell and the model has U = 4, 
V = -0.5, and e/ = -5. 



C. Quantum Monte Carlo Results 

All our simulations were performed for lattices of 4 x 4 
and 6x6 unit cells. The cost of performing these simu- 
lations is approximately the same as simulating an 8 x 8 



We note that we found a PM ground state for one elec- 
tron removed from 1/4 filling, even for large values of U. 
As we will discuss below, when one electron was added to 
1/4 filling, we found a FM state for e/ = —5. Clearly the 
properties of the PAM are asymmetric about 1/4 filling 



in contrast to generally accepted suggestions that it be- 
haves as a half-filled, nearest neighbor hopping Hubbard 
model which displays particle-hole symmetry. 

Figures Q and || show the main results of our simula- 
tions. They plot the energy difference between a polar- 
ized and a singlet ground state as a function of electron 
filling. Figure Q is for the 4x4 lattice, and Fig. Q is 
for the 6x6 lattice. In both figures, t = 1, V = 0.5, 
[7 = 4, and e — —5. For the non-interacting problem, 
these parameters are the "flat band" case illustrated in 
Fig. 0b. 

At 1/4 filling, the two energies are equal to within sta- 
tistical error, indicating the ground state is likely PM. 
This is not what is expected. Our effective one-band 
Hamiltonian has long range hopping. So for a sufficiently 
large value of U, we expect an AF ground-state as was ob- 
served in the one-dimensional DMRG calculations. Cal- 
culating Jrkky for the effective Heisenberg model, we 
found that when i j^ j the magnitude of the exchange in- 
teraction is about two orders of magnitude smaller than 
the magnitude of our statistical error. Thus we should 
not expect to see this state easily. If we used an anti- 
ferromagnetic state for |^t), we would see a AF ground 
state, but the energy of this biased result and the one 
from a PM \iPt) were typically equal within statistical 
error. 
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FIG. 4. Energy difference between the partially polarized 
and singlet ground state energies as a function of electron fill- 
ing. The lattice has 4x4 unit cell and the model has U = 4, 
V = -0.5, and e/ = -5. 

By doping 1/4 filling with one extra election, we found 
that the FM state has a lower energy. Adding one elec- 
tron more produced a PM state. This behavior with the 
doping of one and two electrons is reminiscent of Nagaoka 
ferromagnetism and its instability in the [/ = oo nearest- 
neighbor-hopping Hubbard model. While our model is 
not the same, the proposed mechanism for ferromag- 



netism is very similar: An electron lowers its kinetic en- 
ergy by moving through a ferromagnetic background. In 
our mechanism, the electron lowers its kinetic energy by 
inter-band processes enabling its hopping between two 
occupied /-states anti-aligned with its spin. The need for 
inter-band processes differentiates our mechanism from 
Nagaoka's. 
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FIG. 5. Ground state energy as a function of the value of 
the total spin S^ for several different values of electron fill- 
ing. The lattice has 4x4 unit cell and the model has [7 = 4, 
V = —0.5, and e/ = — 5.The number of electrons varies as 18, 
20, 22, 24, and 26. 

We believe the disappearance of the FM with the sec- 
ond electron is not analogous to the known instability of 
the Nagaoka state but rather is a finite-size effect. It is 
well documented that the shell structure of a finite-sized 
non-interacting problem is reflected in the behavior of the 
filling dependence of the energy in the interacting prob- 
lem ||3l|-|3^. At least to a first approximation, the chem- 
ical potential is constant in a shell and is discontinuous 
between such shells. According to our effective one-band 
model, if we were to add one electron to the 1/4-filled 
case, it would most effectively lower its kinetic energy by 
moving between oppositely aligned /-states by virtually 
hopping through a k = conduction (upper) band state. 
The Pauli principle blocks, or frustrates, this process for 
a second electron of the same spin if finite-size effects 
make the k = (±7r/2, ±7r/2) states of the next shell ener- 
getically inappropriate: The second electron would have 
to enter the k = shell oppositely aligned with the first, 
and the ferromagnetic background of /-states cannot ac- 
commodate the two different alignments. 

The influence of finite-size and inter-band process is 
also seen in Fig. H where we show the energy ratio 
E/\E{Sz = 0)1 as the number of electrons increases from 
18 to 26 electrons, from one closed conduction band shell, 
Ne = 18, through another, N^ — 24, and ending at an 
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open conduction electron shell. At first the polarized 
ground state becomes much lower than the PM state 
and then two become approximately equal at N^ = 24. 
This variation is highlighted in Figs, pk-e which show the 
ground state energy as a function of Sz for different elec- 
tron fillings. For 18 and 26 electrons, the energy per site 
either does not vary within statistical error or increases. 
For 20, 22, and 24 electrons, it has a clear minimum but 
the minimum for 24 electrons is just barely below unity. 
We note that the minimum in energy and the maximum 
in Sz occurs for a half-filled conduction electron shell. 
It is as if a Hund's rule coupling is operative in a non- 
multi-orbital situation. In Fig. |j, Smax is the value of 
Sz for which E/\E(Sz = 0)| is a minimum. If this ratio 
had no clear minimum, then we used Smax = ^(^e — Ng) 
where N^ is the number of lattice sites (the number of 
/-orbitals) and N^ is the number of electrons. 

Nf. — 24 corresponds to 3/8 filling. Instead of being a 
possibly a weak FM state, we argue that it is actually an 
unpolarized state RSDW state; however, while we do not 
clearly see this state for the 4x4 lattice, we clearly saw 
it for the 6x6 lattice when N^ — 56 |£^. There is an 
important and subtle difference between the two lattices: 
Our performance gain in constructing l^r) is based on 
using closed shell states of the non-interacting problem, 
and these states mix information about both the valent 
and conduction bands. For these states, the N^ — 24 
singlet for 4 X 4 is a half-filled shell, but the N^ — 56 
singlet for 6 X 6 is a closed shell. This difference makes it 
much more difficult to see the RSDW state for the 4x4 
lattice. We emphasize that the shell in Fig. is different. 
It is a remnant of a shell in the conduction band. 

We remark that our effective Heisenberg models admits 
a RSDW state. In particular, within the effective Heisen- 
berg model, the ferromagnetic nature of the RKKY ex- 
change interaction depends on the volume enclosed by 
the Fermi surface being small. Increasing the number 
of electrons increases the volume, lengthens the k- vector 
for transfer across the Fermi surface, and decreases the 
strength of the interaction until the Fermi surface be- 
comes unstable because of nesting. Direct calculation of 
the Fourier components of Jj shows that it diverges 
close to nesting. The spatial pattern of this diverging 
Fourier component corresponds the J^ ■ being zero if sites 
i and j are on different sublattices and behaving anti- 
ferromagnetically if the sites are on the same sublattice. 
This is precisely the pattern observed in the spin-spin 
correlation function reported in our previous work, but 
we do not see it for the 4x4 lattice size. Finite sizes 
effects and the limited range of system sizes we can af- 
ford to simulate have prevented us from performing the 
scaling analysis need to establish the RSDW state as one 
of long-range order. 

Both the ferromagnetism and conduction band shell 
structure is evident for the 6x6 case shown in Fig. a. 



Adding one electron to the 1/4 filled case {Ng = 36) pro- 
duces a FM state. Adding another closes the k = (0, 0) 
shell which frustrates the ferromagnetism. Adding more 
electrons to reach N^ = 54 successively populates the 
k = (±7r/3, ±7r/3) and (±27r/3, ±27r/3) conduction band 
shells. At the middle of each shell, the energy difference 
is a maximum and the shell is half filled. 



V. CONCLUSIONS 

Our numerical results indicate that several important 
features of the phase diagram of the one dimensional 
PAM |g,0| are preserved in two dimensions. In both di- 
mensions and for half-filling, the Coulomb interaction in- 
duces an insulating gap, and the system can have AF 
order (AF insulator), or it can remain in a paramagnetic 
state (Kondo insulator) if there is a strong enough hy- 
bridization V between both bands. In this latter case 
there is also a spin gap associated with the energy neces- 
sary to break a Kondo singlet. The AF order originates 
in the Fermi surface nesting at Q = (tt, tt). 

When doped away from half-filling, the system in two 
dimensions becomes paramagnetic. In one dimension, 
however, critical inconmensurate correlations peaked at 
Q = 2k/ develop. This can be understood with the ef- 
fective Heisenberg theory derived in Section IV. B: The 
RKKY interaction J^- ' has a divergent Fourier compo- 
nent at Q = 2k/ because there is nesting for any con- 
centration of electrons. This results in a Luttinger liquid 
with spin-spin correlations which are critical and peaked 
at Q = 2k/. The situation is different in two dimen- 
sions where there is nesting only for fillings 1/2, 3/8, and 
close to 1/4 (small Fermi surface). Our numerical evi- 
dence suggests that the system is paramagnetic between 
3/8 and 1/2 filling. For 3/8 filling, where the nesting 
appears at Q = (0, tt) and (tt, 0), our CPMC results in- 
dicate the presence of a RSDW phase. In this phase, the 
two interpenetrating sublattices are decoupled, and there 
is AF order in each. This spatial ordering can be also un- 
derstood by considering the effective Heisenberg theory 
for the PAM: It is clear that in going from 1/2 to 1/4 fill- 
ing the nearest neighbor RKKY interaction changes its 
sign (AF close to 1/2 and FM close to 1/4). A cancel- 
lation then must occur at some intermediate concentra- 
tion. We can easily see that this intermediate filling is 
3/8 where there is nesting for two different wave vectors: 
Q = (0, tt) and (tt, 0). These two transferred wave vec- 
tors give canceling contributions to the nearest neighbor 
RKKY interactions. For instance, the nearest neigbhor 
in the x-direction will feel an antiferromagnetic interac- 
tion coming from Q = (tt, 0) and a ferromagnetic one of 
the same magnitude coming from Q — (0, tt). On the 
other hand, the next nearest neighbors have an overall 
phase TT for both wave vectors (constructive interference) 
and the effective interaction is therefore AF. This is the 
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origin of the RSDW state obtained for 3/8 filling. 

By decreasing the filling below 3/8, the CPMC results 
indicate that the system becomes paramagnetic again 
down to some critical filling beyond which non-saturated 
ferromagnetism appears. Again this behavior is related 
to the nesting of the Fermi surface. Below 3/8 filling 
there is no nesting down to some small concentration 
of conduction electrons near 1/4 filling where the Fermi 
surface of the conduction electrons can be very well ap- 
proximated by a small sphere and the wave vector 2k/ 
connecting two different points is very close to zero. Un- 
der these conditions, the RKKY interaction J^ ■ diverges 
with a negative value indicating the presence of a ferro- 
magnetic instabilty. 

This mechanism is not the only one giving rise ferro- 
magnetism in the PAM for small concentrations of car- 
riers. In the Mott insulating regime {U < A) where 
the one band effective model Hi band is a valid effective 
low energy theory, the ferromagnetism is related to the 
long range hopping processes involving charge virtually 
tranferring from the lower to the upper band. The rele- 
vance of these processes is related to the flatness of the 
lower band and the dispersive character of the upper one. 
Without considering processes connecting the two differ- 
ent bands {Hiband — -Riband)' ^^^ effective one band 
Hamiltonian is a an extended Hubbard model (the hop- 
ping is not restricted to nearest neighbors). Due to the 
flatness of the lower band, the effective ratio U/t^ is 
very large, and the model is in a Nagaoka-like region even 
for small values of the bare interaction U. The relevant 
low energy scale which determines the magnetic order- 
ing comes from the comparison between pr^j (energy per 
link of the Nagaoka state) and t"j'^/U (energy per hnk of 
the AF state). As this difference is extremely small, any 
other term added to the Hamiltonian, which favors one of 
the two competitive phases, will be relevant. As we have 

(2) 

explained in Section IV. A, H^f^^^^ stabilizes the ferro- 
magnetic phase. If the localized electrons are polarized 
in the same direction, the added carriers can gain energy 
from virtual processes transferring charge between k ~ 
states of the lower and the upper band. 

It is also important to make some additional comments 
about the finite size effects in our numerical results. In 
an infinite system it is necessary to have a finite critical 
concentration of electrons in order to induce the ferro- 
magnetic phase. We argue that This must be the case 
because the gain in kinetic energy of the added parti- 
cles must overcompensate the loss of magnetic energy of 
the localized electrons (which is proportional to the sys- 
tem size). The numerical results show that the system 
becomes ferromagnetic with the addition of only one con- 
duction electron. This is because one electron added to a 
4 X 4 or 6 X 6 system corresponds to a finite concentration. 
Besides, as it is explicitly shown in Figs. and H, there 
are closed shell effects that give a non-monotonic behav- 



ior for the energy difference between the ferromagnetic 
and the paramagnetic state as function of concentration. 
It is clear that for small systems and closed shell con- 
ditions, the non-interacting (paramagnetic) solution will 
be more stable under the introduction of correlations. 
The small systems will generally have larger energy gaps 
between two shells. 

Due to size effects, we must assume that we are deal- 
ing with states of long-range order, and we cannot say 
much about the order of the paramagnetic-ferromagnetic 
transitions. To estrablish long-range order and to deter- 
mine the order of the transitions, it is necessary to scale 
the magnetization with the system size for a fixed filling. 
To do this properly requires larger systems than the ones 
considered in this paper. 
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APPENDIX: A 

To make the derivation of Hspin clearer, we begin by 
rewriting Eq. n^ as 



H = i?(") + i/(i) 



(Al) 



with 



U 



H^'^-T.-.dl'i^. + ^fE-L + TT. 



i i 



ff(i) ^VY. (e^'-'-'dL/.. + e-^'^-'-'/.UkJ 



(A2) 



(A3) 



The next step is the application to iJ of a Schrieffer- 
Wolff-like transformation |3^ to eliminate terms of order 
V and V^ from the Hamiltonian. We do this by expand- 
ing the exponentials of Eq. and imposing the condi- 
tions 



i7(i) + [vi/i,H(")] = o 



(A4) 



;[VFi,[VFi,i/W]] + [M^3,i^(''^] = (A5) 



that define Wi and W^. With these conditions 



(A6) 
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The following expression for Wi satisfies Eq. A4 



Through fourth order in V 



i,k,a (A7) 



with 



and 



7k = 



Ck ~ ef ~U ek - 6/ 



ek-e/ 
With this Wi , the second order term in H becomes 



(A9) 
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J2e^i^~^)-^[t<ff+i^(n^,+n^,)]flU 



■1,3,17 
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where ?Tikk' = 7k + 7k' , the components of the spin op- 
erator Si are 



^i — -^{fi\fi] ^ ftifii) 



Si — fi^fii 



(All) 



and the components of Sg-cr' are the Pauli matrices di- 



f(i) 



,(2) 



vided by two. t]-' and t]-' are defined as 

'■J '■J 
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i^-^j:^''-2s^ 
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Finally Hspm is obtained by means of a second canon- 
ical transformation which eliminates H^'^'> from H 



H — p^ Hp~^ 

^-^ spin — c- -11 1:- , 



By expanding the exponentials of Eq. A13 , we get an 
elimination condition that defines 5: 

ij(2) + [s, H]=0 (A14) 

A S that satisfies this equation is 

?,k,k ,cr,(7 

+ E e*^''-''')--(tij) + i|f )K, - n^,)flh^ (A15) 



iJ,^„ = ij(0) + l[5,ij(2)] + 1^1, ri,i?(2']] 



(A16) 



The second term in Eq. |A16| reduces to 



(A8) i[5,ij(2)] = ]::! ^ e-.(k-k').R, 



X ,^K)o(^-«)o)..L, + ^"i^^^S. ■ S, (A17) 



Ck - ek' 



'kk' 



c/ 



where {n^Q — (i ^^ c^ko-'^kcr)o and (• • • )o means the ex- 
pectation value relative to non-interacting conduction d- 
electrons {V = 0). The first term of Eq. A17| gives the 
RKKY interaction. 



The third term of Eq. |A16| , H'-^'> = i[VFi, [Wi^H^^m 
can be expressed as 



^(4) 



i_ 

4A^2 



i,jMM' 



-i(k-k')-Rij 



[3(5k(5k'TOkk' 



[^k + <5k' - 4mkk' + 2nk(37k + 7k')]ykk'"kk'7k7k' 
4mkk'2/kk' - 77k7k' + 2nk7k7k'(37k + 7k') 
27k7k'™kk']Sj • Sj (A18) 



with 



"^kk' = ^k + (^k' 

a;kk' = ^k + Sk' - 4(7k + 7k') + 2/k(37k + 7k') 

2/kk' = <5k7k' + <5k'7k. (A19) 



APPENDIX: B 



r(2) 



In this appendix we show the derivation of Hi ^^^^^ by 
means of a second order canonical transformation. The 
starting Hamiltonian is: 






(Bl) 



r(0) 



(A13) where ^{bimd ^^ given by Eq. (p5|). From Eqs. ( p6[ ) and 



(1) 



) we can rewrite 7J| j,^^^^ in the following way: 



H 



(1) 

1 band 






The canonical transformation is given by the following 
equation: 



H = e^'He-^^ 



r(i) 



(B3) 



The terms of order H^ j,^^^ are eliminated under the fol- 



lowing condition: 



13 



H 



(1) 

1 band 



r(n) 



[^I'^lfcand] -0 



(B4) 



Satisfying Eq. (gj) is 



Ti ==2 



N\ef\ 



i.k.cr 



£;+(k)-e7-C/ 



Sk) 

^/ (B5) 



where Sk = cos kx + cos A:j,. 

Ill this way we get the second order part of the trans- 
formed Hamiltonian H: 



[9] 

[lo; 
[11 

[12 

[13 

[14 

[15 



1 hand i\j\ ^ \ 



E 

2,J,k,cr 



n" 7?" fo'''' 



k-^f-^ (B6) 



where 



t,,(k) 



1, for i — j 

Y^Sk, for i and j n. n. 

,2 o 

^^tS, , for i and i second and third n. n. 



Here we have considered that the lower band is disper- 
sionless with E~ (k) ~ e/ and e/ = r,^ — r," . 
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